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We have developed a relativistic coupled-cluster theory to incorporate nuclear spin-dependent 
interaction Hamiltonians perturbatively. In this theory, the coupled-cluster operators in the elec- 
tronic sector are defined as tensor operators of rank one and we introduce suitable diagrammatic 
representations. For properties calculations, the electronic part is first calculated and later coupled 
to the nuclear spin part. The method is ideal to calculate parity violating nuclear spin-dependent 
electric dipole transition amplitudes, i51p^C) of heavy atoms. To validate the proposed method the 
£lpNc of the transition 6 ^3^/2 7 ^51/2 in ^^^Cs is calculated for selected MBPT diagrams and 
compared with the results from our theory. 
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I. INTRODUCTION 

Experimental observation of nuclear spin dependent 
(NSD) parity non-conservation (PNC) is a signature of 
nuclear anapole moment (NAM) [1] . It is the most domi- 
nant NSD-PNC effect among three possible sources. The 
other two are: combination of hyperfine and nuclear spin- 
independent PNC; and spin-dependent Z exchange be- 
tween electrons and nucleus. To date, the most precise 
atomic PNC measurement, of atomic Cs, has provided 
the only experimental evidence of NAM 2 . One major 
hurdle to a clear observation of nuclear anapole moment 
is the large nuclear spin-independent (NSI) signal, which 
overwhelms the NSD signature. However, proposed ex- 
periments with single Ba"*" ion [3] could probe PNC in the 
Si/2 — 0^5/2 transition, where the NSI component is zero. 
This could then provide an unambiguous observation of 
NSD-PNC and NAM in particular. The ongoing exper- 
iments with atomic Ytterbium ^ is another possibility, 
the 6s2 ^Sq - 6s5d 3L»2 transition, to observe NSD-PNC 
with minimal mixture from the NSI component. One cru- 
cial input, which is also the source of large uncertainty, 
to extract the value of NAM and nuclear weak charge in 
NSI-PNC is the input from atomic theory calculations. 
In the case of isotope chain measurements, like in Yb, 
the PNC observable is a ratio. Atomic theory contribu- 
tions then cancel and parameters can be extracted with- 
out atomic theory calculations. However, results from 
atomic theory calculations are important in estimating 
the expected value of PNC transition amplitudes and ex- 
tracting NAM does require results from atomic theory 
For these reasons, it is important to employ reliable and 
accurate many-body theory in the atomic theory calcu- 
lations. 

Recently, atomic theory calculations have investigated 
the NSD-PNC of Ba"^ and Ra"*" using the atomic many- 
body perturbation theory (MBPT) 5]. Another recent 
work reports the theoretical estimate of the NSD-PNC 
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observable of Yb and the calculations are based on 
the CI-MBPT method [?]• There is also an earlier work 
on NSD-PNC of Yb using the same method |H] • Very re- 
cently, the CI-MBPT method is used to to calculate the 
NSD-PNC observables of Ba+, Yb+ and Ra+ [9,. Be- 
sides the methods mentioned, earlier works on Ba+ [TU] 
and Yb [11] used configuration interaction (CI). To date, 
coupled-cluster theory (CCT) considered as one of the 
most reliable and accurate many-body theory has not 
been used in the NSD-PNC calculations. The difficulty 
of developing a suitable theory stems from the compli- 
cations of dealing with the nuclear spin dependent in- 
teraction. In a previous paper of ours we reported the 
development of a relativistic coupled-cluster (RCC) the- 
ory to incorporate nuclear spin-dependent interaction in 
a consistent scheme [121 . In this work we provide elabo- 
rate details of the theory and touch upon related subtle 
issues. 

The coupled-cluster (CC) theory [l3l 114] is one of the 
most reliable many-body theory to incorporate electron 
correlation in atomic calculations. In atomic physics, the 
relativistic coupled-cluster (RCC) theory has been used 
extensively in atomic properties calculations, for exam- 
ple, hyperfine structure constants (TSl [H] , electric dipole 
moments |17l I18| , and electromagnetic transition prop- 
erties [T9j[2Q]. In atomic PNC calculations too, RCC is 
the preferred theory and several groups have used it to 
calculate NSI-PNC of atoms pTH25] . However, the cal- 
culations in Ref. [3T] are entirely based on RCC with a 
variation we refer to as perturbed RCC (PRCC), where 
as the calculations in Ref. E5] are based on sum over 
states with CC wave functions. Naturally, the former in- 
corporates electron correlation more precisely than the 
later approach. 

In this work we provide a detailed explanation on the 
approach we have adopted to incorporate nuclear spin- 
dependent interaction Hamiltonian as a perturbation in 
RCC theory. In Section. |llj we give a very short ac- 
count of RCC theory to serve as a quick reference. We 
then give the details of our formulation to represent a 
NSD interaction perturbed CC operator in Section. HI 



PRCC equations of closed-shell and one-valence systems 
is derived in this section. To simplify angular factor cal- 
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culations, the calculations in the electronic sector is sep- 
arated out. Based on this the calculation of NSD-PNC 
observable, in the electronic sector, is described in Sec- 
tion. |IV[ Next section. Section. |Vj describes the coupling 
of results from the electronic sector with the nuclear spin. 
This completes the theoretical development and then, the 
correctness of theory is verified in Section. |VI[ For which, 
we use MBPT results and calculate the NSD-PNC of the 
6 ^5*1/2 — 7 transition in Cs. We then conclude 

and the angular factors of all the diagrams in the lin- 
earized PRCC are given in the Appendix. 



II. BRIEF REVIEW OF RCC 

For A^-electron atoms or ions the Dirac-Coulomb 
Hamiltonian, appropriate to account for relativistic ef- 
fects, is 



N 

E 



(1) 

where cxi and /3 are the Dirac matrices, p is the linear 
momentum, Vn^t) is the nuclear Coulomb potential and 
last term is the electron-electron Coulomb interactions. 
For one- valence systems it satisfies the eigen value equa- 
tion 



(2) 



where j^'i,) and are the atomic state and energy re- 
spectively. In the CC method, the atomic state is ex- 
pressed in terms of T and S, the closed-shell and one- 
valence cluster operators respectively, as 
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1 + ^(0) 



1$,: 
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where is the one- valence Dirac-Fock reference state. 
It is obtained by adding an electron to the closed-shell 
reference state, |<I>i,) — aJ^|$o)- In the coupled-cluster 
singles doubles (CCSD) approximation T^o) = t[°^ +T^°^ 
and are the solutions of the nonlinear coupled equations 



($S|5n|*o) = 0, 

«|iJN|*o)=0, 



(4a) 
(4b) 



where — e "^^ '-/Jnc^' ' is the similarity transformed 
Hamiltonian and the normal order Hamiltonian Hjsi = 
H-{^o\H\<i>o)- The states and \<^ll) are the singly 
and doubly excited determinants, respectively and abc . . . 
{pgr . . .) denote occupied (virtual) orbitals. The details 
of the derivation are given in Ref. m]. Like in T'^°\ the 
one-valence cluster operator 5'^°^ in the CCSD approxi- 
mation is 5"^°^ = s[^^ + S^'' . And these are the solutions 
of 

{^p\Hn + {h7s^°^}\^v) ^ E^''{^P\S[°^\<^^), (5a) 



where = — Eq, is the attachment energy of the 
valence electron. In our previous work \25\ . we provide 
details of the derivation. 



III. PERTURBED CC WAVE FUNCTION 

Time independent perturbation theory is the standard 
procedure to incorporate external perturbations or ad- 
ditional interactions in atomic many-body calculations. 
However, a basic requirement of perturbative calculations 
is a complete set of intermediate atomic states, which 
is non-trivial to generate. The perturbed CC method 
[Ml [33 [23 J on the other hand, implicitly accounts for 
all the possible intermediate states. In this work, we 
consider the perturbation as the nuclear spin-dependent 
parity non-conserving (PNC) interaction 
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PNC 
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where ^'^y is the weak nuclear moment of the nucleus 
and p-^ir) is the nuclear density. The weak nuclear mo- 
ment is expressed in terms of the neutron and proton 
numbers = 2(ZCip -I- NCin), where Cip and Ci„ 
are respectively the vector electron and the axial vector 
nucleon coupling coefficients. There are two complica- 
tions arising from the nuclear spin operator / in H^^. 
First, the cluster operators in the electron space are rank 
one operators, and second, the atomic states in the one- 
valence sector are eigenstates of total angular momentum 
F = I -I- J. Both of these are relatively simple to incor- 
porate at lower order MBPT calculations, however, im- 
plementing these in non-perturbative theory like RCC is 
nontrivial. Accordingly, the method we have developed 
and implemented in the current work are very different 
from our previous works [241 12 7| . The other simplifying 
feature of the previous works is, the methods developed 
were for closed-shell systems where the total electronic 
angular momentum J = 0. 

With the PNC interaction, the total atomic Hamilto- 
nian is 



Ha = H 



DC 
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where A is the perturbation parameter. As Hp^^ mixes 
atomic states of opposite parities, the eigenvalue equation 
is modified to 



HA\'fv)=EJ-^, 



(8) 



Note that the first order energy correction E^ = 
(^'oliJplgl^'o) = as i?pN§ is an odd parity interaction 
Hamiltonian. This is taken into account while writing 
the above eigenvalue equation. Here l^^v) are the mixed 
parity atomic states and to first order in perturbation 
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E„ — Et 



(9) 
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If \^v) and l^'io) are atomic states of same parity, then 
the iypNc induced electric dipole {El ) transition ampH- 
tude is 



El 



NSD 
PNC 



= (*,„||D||*,, 



(10) 



where D is the dipole operator. The perturbation ex- 
pression in Eq. (|9| require a complete set of intermediate 
atomic states This, as mentioned at the beginning 

of the section, is non-trivial to obtain in atomic many- 
body calculations. Summation over intermediate states is 
circumvented when E'lp^p is calculated with CC atomic 
states. For this define a new set of cluster operators T^^^ , 
which unlike T^*^^ connects the reference state to oppo- 
site parity states. This is the result of incorporating one 
order of Hp^Q and for this reason we refer to T^^^ as the 
perturbed cluster operators. Although hyperfine states 
are natural to _ffp^§, cluster operator T^^^ is defined to 
operate only in the electronic space and is a rank one 
operator. So, the mixed parity atomic state in RCC is 



where \ '^v) is the opposite parity component arising from 



the iJpj^Q. The first order wave function is then, 



m = 



1 + n 



Although, notation wise il^^J seem first order, it is sec- 
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ond order in perturbation: 



Coulomb interaction and Hp^^ 
sisting of same parity states 



one order each in residual 
For model space con- 
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Here, and P and Q are the projection operators of the 
model and complementary spaces, respectively, Hq is the 
Dirac-Fock Hamiltonian, Vfcs is the residual Coulomb in- 
teraction and the zeroth order mixed parity wave opera- 
tor 
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The scalar product with the nuclear spin I in the expo- 
nent restores T^^^ to the correct form of the wave oper- 
ator. At this point it is convenient to separate out the 
electronic part of the interaction Hamiltonian 
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elec 
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which operates only in the electronic space, so that 
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The remaining part of this section describe how to arrive 
at a consistent representation of T*^^-' and extension of 
the method to one-valence systems, where S*-^^ are the 
perturbed cluster operators. 



Using fl'^J'l and fl^l, we may compute SlpNC to third 
order in perturbation, the details of which are discussed 
in Ref. ^5j. 



B. Perturbed cluster operator representation 

The expression of fl^'^l and associated selection rules 
are what we need to arrive at a consistent description 
of T(i). In Fig. [ija) we show an MBPT wave function 
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A. MBPT wavefunctions 

To define the multipole and parity selection rules of 
T(i) , we examine the second order MBPT wave function 
with -ffpNc ^-s one of the perturbations. From the gen- 
eralized Bloch equation [28], the total wave operator is 

m 
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in) 
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Here fl^ J, has n orders of residual Coulomb interaction. 



'u,0 

where as 17, 
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v,l 



has one order of Hp^^ and n orders of 



residual Coulomb interaction. Following which the mixed 
parity state 



I*, 



I*,: 



I*, 



FIG. 1. Second order MBPT diagram where interaction 
Hamiltonian i^pNc operates on the particle state p. (a) Dia- 
grammatic representation of Hp§j^ in the electronic and nu- 
clear spin space, (b) Diagrammatic representation within the 
electronic space only where the PNC interaction Hamiltonian 
is H^i^^ . Line terminated with square represents the PNC 
interaction. 



diagram of fl^'^l arising from the second term in Eq. ( 18 1 
The algebraic expression of the diagram is 
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where are the orbital energies. The states | • • • /) rep- 
resent uncoupled atomic states comprising of electronic 
and nuclear states. 

On closer examination, the component of iJpNc which 
operates in the nuclear subspace is diagonal. So, for fur- 
ther calculations, we separate out the matrix elements 



FIG. 2. Angular reduction of the second order MBPT dia- 
gram. The remnant angular momentum diagram on the right 
is representative of the multipole structure of PRCC double 
excitation cluster operators Tj'^' and Sj^'. 



FIG. 3. Diagrammatic representation of the single and double 
excitation perturbed cluster operators in closed shell and one- 
valence sectors. The extra line in the Tg^' and Sj^' is to 
indicate the multipole structure of the operators. 



in the electronic subspace and combine with the nuclear 
part at the end of the calculations. For the MBPT di- 
agram considered earlier Fig. [ija), the electronic part 
is as shown in Fig. [ijb). The corresponding algebraic 
expression is 



{w\oLp^{r)\p) {pq\V\va) 
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Details of the angular reduction when electronic state is 
coupled with the nuclear state are discussed in the later 
sections of the paper. Like wise, the perturbed cluster 
operators T^^\ as defined earlier operate only in the elec- 
tronic. For this consider the angular part of the matrix 
elements in Eq. (20), the diagrammatic representation 
is shown in Fig. [2] The angular diagrams are based 
on the conventions used in Lindgren and Morrison |28j 
and the same is followed in the remaining of the paper 
while referring to angular momentum diagrams. In the 
figure, the angular diagram on the right hand side in- 
dicates the multipole structure of the t!^^ and like the 
electron-electron Coulomb interaction 



^2^^ = EE^a7a,fc2){Q(fi)C,,(f2)}\ 

abpq l,k2 



(21) 



where Ci are c-tensor operators and {• • • }^ indicates the 
two c-tensor operators couple to a rank one tensor op- 
erator. Following general rules of coupling tensor op- 
erators, the rank of the tensor operators must satisfy 
the triangular conditions \jw — jv\ ^ I ^ [jiv + jv), 
\ja - jq\ fc2 < {ja + iq) and |Z - /caK 1 + k2). 
Effectively, T*^^-' is a rank one operator in the electronic 
subspace and forms a scalar operator after coupling with 
I. The other important difference from the T^'^^ is the 
parity selection rule at the vertices. The combined pari- 
ties at the vertices are opposite (—1)'™+'" = _(_i)'<i+'g 
. After a similar analysis, the singles operator is 



.(1) 



(22) 



Based on the multipole structures, the perturbed cluster 
operators are diagrammatically represented as shown in 
Fig. [s] For the doubles Tj^'', to indicate the multipole 
structure, an additional line is added to the interaction 
line. 



C. Closed-shell systems 

The mixed parity states, eigenstates of i?A, in the per- 
turbed relativistic coupled-cluster (PRCC) is as given in 
Eq. (11). However, as i?pNc is considered to first order 



only, it is sufficient to consider up to linear terms in the 
expansion of e^'^^ The PRCC atomic state is then 
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i + at(i) - I 
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The cluster operator T^^^ as described earlier, incorpo- 
rates H^j^iP to first order and residual Coulomb interac- 
tion to all order. The closed-shell equivalent of Eq. (|8|, 

the mixed parity eigenvalue equation, is 
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From the definition of the coupled-cluster mixed parity 
state in Eq. ( 23 ) , the eigen value equation in terms of 



the coupled-cluster wavefunction is 
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AT(i) • I 
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In terms of Hjss, defined earlier, the eigenvalue equation 
is simplified to 



AH 



NSD\ 



= AEoe 
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where AEq = Eq — {^o\H\^o) is the closed-shell correla- 
tion energy. Retaining terms which are first-order in A, 
we get 



iJNT(^) • I - 



HNSD T 
clcc ■ 



e' l^o) 



= A£;oT(i).Ie^*'"|$o; 

(27) 

From here on, for simplicity, we drop the nuclear spin 
operator I and work only in the electronic space. At 
a later stage of the calculations the electronic part is 
coupled with the nuclear spin. Operating with e~"^' ' and 
projecting on singly and doubly excited states and 
(<i>^'|, respectively, we get the CC equations for singles 
and doubles perturbed cluster amplitudes as 

($S|{H^(i)}|$o) - -($S|HSe?l$o) 
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FIG. 4. Diagrams which contribute to the linearized per- 
turbed coupled-cluster equations of the singly excited cluster 
operator Tj^^\ 
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FIG. 5. Diagrams which contribute to the closed-shell lin- 
earized perturbed coupled-cluster equations of the doubly ex- 



cited cluster operator T, 



(1) 



Where, H^SD ^ e-T(o)H. 



'e" , is the similarity trans- 



NSD T<°' 
dec ^ ' 

formed PNC interaction Hamiltonian in the electronic 
subspace. The perturbed cluster operators are solutions 
of these coupled linear equations. However, the equations 
are nonlinear in the unperturbed cluster operators T^^\ 
The advantage of separating the cluster operators into 
two categories T^"^ and T^^' is, the two sets of equations 
can be solved sequentially. As the T^^^ has no dependence 
on T*^^\ the T^^^^ equations are solved first and then the 
T*^^-' equations are solved. An approximate form of Eq. 
(28a I and (28b), but which contains all the important 



many-body effects, are the linearized cluster equations. 
This is obtained by considering 
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The diagrams in the singles and doubles equation arising 



from HnT^^") are shown in F 



ig.HandQ 



respectively. 



The other form of Eq. ( 29 ) is to write the equation in 



terms of specific cluster amplitudes. The singles equation 
is then, 



hr, 
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"^-^gbcqaT^^, (30) 
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where gijki and hij, are the matrix elements {ij\V\kl) 
and (zIH^^jPIj), respectively, and and t^^^ are the un- 



perturbed single and double excitation RCC amplitudes, 
respectively. For compact notation, we have defined 
gijki = gijki - gijik- Similarly, tjfj and t^-^ are the anti- 
symmetized unperturbed and perturbed CC amplitudes 



respectively. The equation of the double excitation per- 
turbed cluster amplitudes is 
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+ Y dpqrbTa - Y 9cqabT^ + ^ 

r c rc 
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(31) 
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where ( ) represents terms similar to those within 
parenthesis but with the combined permutations p O g 
and a b. These equations are similar to the unper- 
turbed cluster equations in Ref. |28) . There is, however, 
a major difference from the all-order equations of un- 
perturbed cluster operators used in other works [211 130] 
which use antisymmetized cluster operators. Here, all 
order is the same as the linearized coupled-cluster theory 
and in the antisymmetized representation = —t^^. 
For our work, which is based on diagrammatic evalua- 
tion of angular factors, the representation without anti- 
symmetrization is preferable as it follows directly from 
the diagrams. To solve Eq. (30) and (31), the expres- 



sions of are separated into multipole components and 
the angular factor of each terms in the two equations are 
given in the Appendix. 



D. One- valence systems 



From Eq. ([24| , we may write the perturbed eigenvalue 
equation of the one-valence system as 



DC 



AHl^.^ . I) Ivp, 



(32) 
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As defined earlier, Ey is the energy of the one-valence 
system. Similar to the closed-shell case, the perturbed 
wavefunction in CC theory is 



1-1- AT 
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1 + S^"^ + AS 
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(33) 

where S'^-* is the perturbed CC operator of the one- 
valence part. The diagrammatic representation of the va- 
lence single and double perturbed cluster operators sj^"* 
and S'^\ respectively, are shown in Fig. [3] These are 
topologically equivalent to the diagrams of the closed- 
shell operators T^^"* and t'^\ however, with one of the 



core lines rotated. The Eq. (33 1 in terms of the PRCC 
wavefunction is 



l-f-AT(i) - I 



1 + 5(0) 
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To derive the S'^^ equations, like in the closed-shell case, 
project the above equation on e^'^ and retain the terms 
linear in A. For further simplification, use normal-ordered 
form of the Hamiltonian, which for the one-valence sys- 
tem is = H — {^y\H\^y) . After these sequence of 
operations, and retaining only the electronic part like in 
the closed-shell case, the eigen value equation is modified 
to 



iJNS(i) + FnT(i)(1 + S) + HNSD(i + s) 



A£;,S(i) + AS„t(i)(1 + 5)1 |$.),(35) 



where AE^ — Ey — is the correlation en- 



ergy of the one- valence system. Projecting Eq. (35 1 with 



the excited determinants ($^1 and (<&?!a|, we get the per- 
turbed CC equations of the singles and doubles respec- 
tively, in the form 
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Where, we have used the relations 

($p|T(i) 1$,) - 0, and ($p|T(i)5|$,) = 0, (37) 

as T*^^-', being the cluster operator of closed-shell sector, 
does not contribute to the CC equation of Sl^"* and Sj^'. 
In terms of specific components, as in Eq. ( 30 1 , the single 



excitation perturbed cluster amplitude equation is 
{ey+AEy-ep)TP = hp^ + hpqtl - ^ hb^tl 
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This is Eq. (30 1 with two important modifications. First, 



the valence orbital v replaces the orbital a and second, 
the single particle energy of v is modified to include the 
correlation energy. Following similar modifications, the 
equation of the double excitation cluster amplitudes is 



{ey+AEy+€b - €p- eq)T^b = E ^P^'^vI ^ E 
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Angular factors obtained from the closed-shell sector, 
with replacement of a by v, can be used to rewrite the 
cluster equations in terms specific multipole components 
from these equations. With these modifications, the same 
methods used to solve the T^^^ are used to calculate the 
5^^^ cluster amplitudes. The fore going derivations shows 
that PRCC equations of one-valence sector are not very 
different from the closed-shell equations. But the same 
rationale does not apply to two- valence PRCC equations. 
For two-valence systems, even at the level of RCC, re- 
quire due considerations on the nature of model space. 
Properties calculations from the RCC wave functions is 
another level of sophistication, these and related issues 
on RCC calculations of two-valence systems is explored 
in one of our recent works |31j . 



IV. Ei: 



CALCULATION 



In this section, as prelude to the calculation of -El^ftP 
from the PRCC states, the details of i^lp^^ calculation 
with MBPT is discussed at the beginning. The MBPT 
calculation, however, is in the hyperfine states and is 
equivalent to the PRCC expressions of El^^^^ after cou- 
pling with the nuclear states. 



From MBPT wavefunction 



To calculate the Hp^^ induced electric dipole tran- 
sition amplitude ElJ^^ with MBPT, we use the wave- 



function in Eq. (16). For the purpose of using MBPT 
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calculations as the basics of analyzing the PRCC, it is 
sufficient to consider the total wave operator 

n = n^^l + n^^l (40) 

The PNC induced El transition amplitude is then 
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The expressions arising from these terms contain all the 
symmetry information to describe the properties of T^^\ 
Consider terms linear in A, 



(42) 



This is the MBPT expression of i?lp^c, which has one 
order each of Hp^^ and residual Coulomb interactions. 

The diagrams which arises from Dfi^^J and ri^^Q^Dfi^^J 
are shown in Fig. |6| A similar set of diagrams, which 
t — t 

arise from D and Dl]|,^^, are obtained from 

interchanging Hp^^ and D vertices in the diagrams. 
Besides providing insights on the nature of perturbed 
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FIG. 6. -^1^1^^*^ diagrams which arise from the MBPT terms 
Dn[°> and ni^^ Dn^"}. These are the first and third terms 



of the Eq. ( |42[ ). In the diagrams interaction lines terminated 
with circle and rectangle represent D and -ff^^'"', respectively. 



amplitude is 



Tt 



1 + AT(i)-I 1 + 5 + AS(i)-I 



D 



xe 



i + at(i) - I 



1 + 5* + AS 



(1) 



ll'i>.)- (43) 



Retain terms linear in A as the remaining are zero from 
parity selection rule. Considering only the electronic 
component, like in the calculation of the cluster ampli- 
tudes T(i) and S(i\ define as the i/^^SD induced 
El amplitude in the electronic space. From the commu- 
tation relation of the cluster-operators, the expression is 
reduced to 



piNSD 



(*^||D 



T(i)+sW 



T(i) + S(i) + T^^'>s\ DS'll*^), (44) 



where D — e^^De^, is the dressed electric dipole oper- 
ator. It is evident that D is a non-terminating series of 
the closed-shell cluster operators. It is non-trivial to in- 
corporate T to all orders in numerical computations. For 
this reason D approximated as 



D « D + Dr(o) + t(")'d + r(")^DT(°) 



(45) 



This captures all the important contributions arising 
from the core-polarization and pair-correlation effects. 
Terms not included in this approximation are terms 
which are third and higher order in T^^\ The expres- 
sion used in our calculations is then 



piNSD 



^^(o)tj3^(i) ^ T(i)tj3^(o) ^ j3s(i) + s(i)^D 
+^(o)t]3s(i)+sW^D5(")||$„). (46) 



The Eq. ( |46[ ), as mentioned earlier, includes term up 
to second order in cluster amplitudes. From our pre- 
vious study of properties calculations 25J, we conclude 
that the contributions from the higher order are negligi- 
ble. Selected diagrams from the leading-order and next 
to leading order terms of Eq. ( 46 ) are shown in the Fig. 

El 



cluster operators, detailed MBPT calculations serve an- 
other valuable purpose. These serve as reference calcula- 
tions to validate the PRCC implementations. 



B. From PRCC wavefunction 

Like in MBPT, from the expression of the PRCC wave 
function in Eq. ( 33 I 



the Hp^^ induced El transition 



V. COUPLING WITH NUCLEAR SPIN 

The PRCC diagrams of -Bl^SD Fig. [7] are, as ex- 
plained earlier, defined in the electronic subspace. To 
calculate iJlp^^, the nuclear spin part of the operator is 
coupled with iSl^^^p. The i?lpNc diagrams must then in- 
clude a vertex which operates on the nuclear spin states. 
This is evident from the example A4BPT diagram in Fig. 
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(c) 



(d) 



(E 



(E 



(g) 



(i) 



(j) 



FIG. 7. Few of the leading order PRCC diagrams which con- 
tribute to the Elli^^ of one-valence atoms. 



[TJa). Although complicated, the coupling with the nu- 
clear part of the operator require only angular integra- 
tion, where as the evaluation in the electronic space in- 
volves radial integrals. In our calculations, we use di- 
agrammatic to carry out the angular integrals and as 
mentioned earlier, we follow the conventions of Lindgren 
and Morrison [25] . 

To elaborate on the coupling schemes, we discuss spe- 
cific cases involving single and double PRCC operators. 
For the former, we select the topologically simplest, and 
for the later, one of the more complicated one. So that, 
these represent the range of topology and complexity of 
the diagrams. In each of these the main objective is to 
reduce the electronic part to a form common to a set of 
diagrams. And then, combine with the nuclear part. 




Ju 



jp h 

F F T 

jp jv h 

I I F„ 



11! 



-o ki 



FIG. 8. Angular factor calculation of the -Blp^^ diagram 
shown in Fig. 71a) and it arises from the term DS[^\ The 
portion within the rectangle in dashed line indicate the elec- 
tronic part. Remaining portion is from the nuclear spin and 
hyperfine coupling. 



A. Singles PRCC operator 



Fig. [7|^a). It represents the Dirac-Fock contribution and 
is perhaps the simplest diagram. But it is the most dom- 
inant and naturally the most important. It is one of the 
diagrams arising from DSj^"* in the PRCC calculations. 
The angular diagram for the electronic part of this dia- 
gram has same topology as in Fig. ^a). However, the 
orbitals lines are replaced by angular momentum lines 
and arrows on appropriate lines to represent phase fac- 
tors. The angular momentum diagram, which includes 
coupling with the nuclear spin, is shown in Fig. [8] The 
electronic part is the portion enclosed with the dashed 
line rectangle. 

In the diagram, and jw are the total angular mo- 
mentum of the single particle states of the initial and 
final one-valance atomic states, respectively. The line ki 
represents a rank one multipole line and denotes the an- 
gular part of the dipole operator. The multipole line ^2 
is the all important representation of Hp^^, it is rank 
one multipole in electronic part as well as the nuclear 
part. In the electronic part, it represents the angular 
part of operator a and produces a transition from jy to 
jp. Here, jp is the total angular momentum of the in- 
termediate single particle state. On the other hand, in 
the nuclear sector, k2 represents the operator I. Since it 
is diagonal operator of the nuclear spin state |/), the ^2 
line does not change the nuclear spin. 

The last step in the angular momentum diagram repre- 
sentation is to coupled the electron and nuclear momenta 
to form the total angular momentum of the hyperfine 
states. The initial and final hyperfine states are \Fy) and 
respectively. The lines marked as Fy and repre- 
sent the angular momenta of the hyperfine states in the 
diagram. Using Wigner-Eckert theorem, we can write 
algebraic equivalent of the diagram as 



TO, 



I\Fymy) - (-1)^ 



F, 



1 Fy 
q m„ 



■{FjD,s\\Fy) 



(47) 



where D^s = DSj^^ • I is a rank one operator, to^ are 
the hyperfine magnetic quantum numbers and q is the 
component of -Dcff- The phase factor and 3j-symbol in 
the above expression are the free lines in the right hand 
side of the diagrammatic equation in Fig. [Sj Remain- 
ing expression on the right hand side, phase factor and 
6j-symbols, is the to independent angular component of 
{Fw\\F>cs\\Fy) . 



B. Doubles PRCC operator 



(1) 



(1) 



For the coupling with nuclear spin involving either T\ 

or Si^, consider a more complicated diagram as shown in 
Fig. [9]ja). It is the exchange of diagram in Fig. [7][j) and 

arises from the term Tj"'' £'82^'' in the PRCC expression 
of -E^J^P. The angular momentum diagram is shown in 
Fig. ^b) and note that the angular mome ntu m represen- 



To define the coupling of singles diagrams, either 
or S^i \ with nuclear spin part consider the diagram in tation of Sj^'' is as described in Section 



III 



The labels 
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FIG. 9. Angular momentum diagram representation of a 
PRCC diagram, (a) Exchange of the dia gram shown 



Fig. 0;,] 

gram. 



(b) The corresponding angular momentum dia- 



of the lines are different from the one in Fig. [2j however, 
these are dummy labels and same selection rules apply 
to h and I2, and fc2 represents a rank one multipole line. 
It has the same role in the electronic sector as fc2 of the 
diagram in Fig. [8j 




10. Angular momentum diagram of SlpNc in the hyper- 



FIG. 

fine states involving double excitation cluster operator ^2 



(o)t. 



The diagram arises from the term 5*2"'' DS'2^' and the por- 
tion within the rectangle in dashed lines indicate the portion 
arising from the electronic sector. 



To demonstrate the non-trivial angular integration in 
the calculations with hyperfine atomic states, the angu- 



lar momentum diagram of (fu,mu;|T'2°'' 03)2^' ■ I|/^m.u), 
of the diagram in Fig. |9l is shown in Fig. jTo] The por- 
tion of the diagram within the rectangle in dashed-line 
is the angular momentum part of the electronic sector 
and except for the topological rotation of the (li, I2, ^2 
) vertex, it is identical to the diagram in Fig. |9jb). The 
evaluation of the angular integral of the electronic sector, 
for the example considered, is shown in Fig. Impor- 
tant point to be observed is the structure of the free part 
in the right hand side of the diagrammatic equation in 
Fig. [TT] Although the multipole lines oi D and a, ki and 
/c2 respectively, are coupled to an effective multipole I4, 




fci = E,„,,.fe.;4](-i)'»+'-' 



f jr h iq \ 

\ h jb hi J 



jw h jp 1 / '4 k2 ki 
jb ki jr j \ h h h 



h h li 

jw jv jp 




FIG. 11. Angular factor reduction in the electronic sector. 
The angular momentum diagram on the left hand side is the 
portion within rectangle in Fig. |10| On the right hand side, 
the angular momentum diagram with free lines represents the 
effective part coupled to the nuclear spin. 



the ki and ^2 are present as free lines. The effective mul- 
tipole I4 operates on jy and transforms it to jy^ . In terms 
of Wigner-Eckert theorem, the diagram is equivalent to 



(j.liTf) DS^l II (48) 



jw ^4 jv 

-ruyj q ruy 



where {. . represents coupling of rank one tensor oper- 
ators D and S*^^-' to an operator of rank I4. This coupling 
is a structure common to any PRCC term of iJl^^^p. 
That is, for any term, the angular integral in the elec- 
tronic sector is reducible to a form where the free lines 



is similar to the one on the right hand side of Fig. 11 



From the triangular condition, I4 = 0,1,2 are the al- 




FIG. 12. Angular factor reduction of the coupling between 
the electronic sector with the nuclear spin of a diagram aris- 
ing from S^''^2DS2^\ The electronic sector of the diagram is 
shown in Fig. |11| 



lowed values, however, what values of 1 4 contribute de- 
pends on jy and jy,. For example, ^4 = 0, 1 contribute in 
the PNC 6 7 ^S'1/2 transition of atomic Cs ^, 
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where as only Z4 = 2 contributes to the proposed PNC 
6 — )■ 5 transition in Ba+ [3]. 

The form of the free hnes in the Fig. [TTJrequire one due 
consideration while combining with the nuclear part. It is 
the multipole line fc2, inherited from T*-^' or S^^-*, which 
couples with the operator I to form a scalar operator. 
The diagrammatic representation is shown on the left 
side of Fig. [12] After evaluation, it reduces to a 9j- 
symbol and free line part. From Wigner-Eckert theorem, 
the matrix element in the hyperfine states is 



F 1 F 

-mm q rriy 



where 



off 



(1) 



(49) 



(50) 



is the effective dipole operator in the hyperfine states. It 
is of the same form as the effective operator in Eq. (47 1. 
However, the coupling sequence in Eq. (49 1 is general 
and applies to all the terms in the iJl^S'^ 



VI. VALIDATION OF PRCC 

The RCC method described so far involves intricate 
but tractable angular momentum coupling and most of 
the calculations are in the configuration space of reduced 
dimension, namely of the electrons. Considering the 
complexity of the method, it is desirable to validate the 
method with few selected terms before a full scale imple- 
mentation. Here, we present a method of validation by 
comparing with dominant third order MBPT diagrams. 
This is possible as we solve the PRCC equations itera- 
tively with the first order MBPT wave functions as the 
initial guess. In particular, we evaluate the £lp^^ ^^i*^ 
transition 6s — 7s in ^^'^Cs. 



TABLE I. Validation of singles PRCC amplitudes. The values 
listed are in units of ieaofj,'^. Numbers in the square bracket 
represent the power of 10. 



Orbital {p,q) 
Gpi/2 

7Pl/2 
8pi/2 
dPl/2 

lOpi/2 

llPl/2 



MBPT 
3.80237[-15] 
-4.42335hl6] 
1.08456[-16] 
1.88997[-16] 
8.66134[-17] 
2.41491[-18] 



PRCC 
3.80237[-15] 
-4.42335hl6] 
1.08456[-16] 
1.88997[-16] 
8.66134[-17] 
2.41491[-18] 



(1) 



A. Single excitation operator S{ 



To check the PRCC equation and angular momentum 
coupling in the calculations of s[^\ consider the term 



S[^^ D in 

grams, which are equivalent to the SY' D at the first 
iteration of the PRCC equation are shown in Fig. [13] In 
the figure, the portion of the diagrams within the rect- 
angle of dashed-line arises from the second order wave 

operator ^^Ji define in Eq. (17). As for the solu- 
tions of linearized PRCC equations, the initial guess is 
the first order MBPT wave function ^'^Ji- Single excita- 
tion cluster amplitudes r are then calculated iteratively 
from Eq. (30). At the first iteration of the linearized 



pNSD 



Two of the third order MBPT dia- 



PRCC equation, the diagram within the dashed rectan- 
gle in Fig. [13]^ a) is equal to the one valence version of 
the diagram in Fig. ^c). And, the diagram Fig. [l^b) 
is the exchange counterpart. These MBPT diagrams are 
the dominant ones after the DF and other theories like 
the PRCC must be able to reproduce matching results. 

For comparison, Slp^^ contribution from the two di- 
agrams m Fig. ]4]from MBPT and resufis fr om the equiv- 
alent PRCC calculations are listed in Table. ]l] The spe- 
cific orbital wise contributions, for better comparison, of 
the dominant contributions are listed in the table. It 
is evident that there is excellent agreement between the 
MBPT and PRCC results. 



7si/2 



6si/2 




' 6si/2 



(a) (b) 

FIG. 13. MBPT SlpNc diagrams used in the validation of 
the the singles PRCC amplitudes S[^K 



7si/2 



6si 



/2 



(a) 



P 



□ 

6si/2 



7si/2 



(b) 



FIG. 14. MBPT Blp^c diagrams used in the validation of 
the doubles PRCC amplitudes S2^\ 
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(1) 



B. Double excitation operator S, 



Like in the case of s[^\ consider two MBPT diagrams 

which are equivalent to the 82^"^ diagram after the first 
iteration and are shown in Fig. |14[ In the figure, the 
portion of the diagrams within the rectangle in dashed- 
lines are equivalent to the 5*2^^ diagram in Fig. but 
adapted for one- valence systems. However, there is a ma- 
jor difference from the s[^K As -ffpNc a single particle 
operator, the double excitation cluster operator at one 
order of TJ^^g is zero and the initial guess value is set to 
zero. To provide a wider test sample, for the double ex- 
citation the MBPT diagrams considered are equivalent 
to the term -DS'j^^' in PRCC. Another variation is, the 

two diagrams considered in sl^"^ arise from topologically 
differen t di agrams, in the present case, the two diagrams 



14 



arise from the same Si^'' but different cluster 



in Fig. 
amplitudes. 



TABLE II. Validation of doubles PRCC amplitudes. The val- 
ues listed are in units of ieaofi'^. Numbers in the square 
bracket represent the power of 10. 



Orbital (p, 



MBPT 



PRCC 



6pi/2 
7pi/2 
8pi/2 
9pi/2 
lOpi/2 
llpi/2 



6.53322[ 
2.35102[ 
3.20651[ 
2.02539[ 
8.67435[ 
3.02443 



-17] 
-17] 
-16] 
-15] 
-16] 
-17] 



6.53322[-17] 
2.35102[-17] 
3.20651[-16] 
2.02539[-15] 
8.67435[-16] 
3. 02443 [-17] 



tum couplings at all stages of the calculations. The pro- 
posed scheme, on the other hand, introduces the nuclear 
spin coupling is with a simplified effective operator in 
the electronic sector. The validity of the representation 
is explicitly tested and verified with selected diagrams. 
Although limited in number, the diagrams and terms se- 
lected in the sample test calculations are varied enough 
to account for complex as well as subtle issues related 
to the method. Based on the results presented, we con- 
clude that the method works and in future publications 
we shall report results of sophisticated and large scale 
calculations using PRCC. 
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Appendix: Angular factors of PRCC equation 

Here, we give the linearized PRCC equations of the 
closed-shell sector, Eqs. (30) and (31 1, after angular in- 
tegration. 



The one-valence cluster amplitude equations 
can be obtained after suitable modifications. 



1. Single excitation cluster operator 



The results from the MBPT and PRCC calculations 
are listed in Table. [H] Here too, the results from the two 
calculations are in excellent agreement. Perhaps, it must 
be mentioned that, the angular momentum factor cal- 
culations of the MBPT and PRCC diagrams are done in 
very different steps. Electronic part of the portion within 
the rectangle in dashed-line in Fig 



reduce to the representation of 5*2^' as discussed in Sec- 
tion. |III B[ MBPT diagrams, on the other hand, may be 
evaluated without the need to associate with an effective 
operator of specific form. 



VII. CONCLUSIONS 

The relativistic coupled-cluster method we have devel- 
oped to incorporate a nuclear spin-dependent perturba- 
tion is an apt one to calculate nuclear spin-dependent 
parity non-conservation in atoms and ions. The repre- 
sentation of the cluster amplitude within the electronic 
sector as a tensor operator of rank one and coupling with 
the nuclear spin part at a later stage of property calcu- 
lation leads to simplification of the calculations. Oth- 
erwise, the entire calculation must be done with hyper- 
fine states, which involves complicated angular momen- 



The angular reduction of each diagram is such that 
the free parts are reduced to the form in Fig. [sja). The 
free part is common to all the diagrams and is avoided in 
the computational implementation. So, solutions of the 
PRCC equations are the cluster amplitudes in reduced 
matrix form. For T^^^^ and s[^\ the rank of the operator, 
14|are evaluated to g^j^^j hence the free part, is one. For convenience, the 



representation of Tj^"* is redefined as 



.(1) 



2222<tihJ2){CiAh)c,,{h)V. (A.i) 



abpq Zi ,^2 



Like the PRCC operators, we use and t^l to represent 
the unperturbed single and double excitation amplitudes. 



respectively. The tensor operator structure of Tj*^^ is 



abpq k 



(A.2) 



For the Slater integrals, the reduced matrix element is 

Xkiabcd) = {-!)'' {Ka\\C''\\K,){fib\\C%d)Rkiabcd), 

(A.3) 

where, Rk{abcd) is the radial part of the Slater integral. 
With these definitions, the Eq. ( 30 ) is written in terms of 



reduced matrix elements and appropriate angular factors. 
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V(^k, + 1) 



(^(fcl, fc2) 

V(2fci + 1) 



E E^.a(^l''2) 

bqrk2 hh 



v/(2A:2 + 1) 



E E^'^^l^cga) 



^(fc2,^2) 

V(2fc2 + 1) 



(A.4) 



2. Double excitation cluster operator 



For the double excitation PRCC operator 
common free part of the angular factors in Eq. (31) is 



the 



as shown in Fig. [sjjb). In the angular factors, the multi- 
poles fci, fc2, li and I2 have the same interpretations as in 
the singles. Multipoles and 1^ are, however, arise from 
coupling {ja,jp) and {jb,jq), respectively. The double 
excitation cluster equation obtained from the projection 
on ($^^| along with appropriate angular factors is 



E(2;i + l)(-l)^"+^^+^-^+'^ I ^1 g > I KrCbih) - J2i^h + 1) 



X (-1) 



Ja+Jp+h / kl ja jc 

jp h h 



- E(2^i + 1)(-1)-+-+^^^-'^ { I g I } X,(c,a6)r,^ + E (^(-D^^^^'^"'^ 
xXfc,(pcar)T;,''(;i,Z2)-EE(2^2 + l)(-iy^-^'+'=^+'-^<| Jc ^4 3q \ Xi,{vcar)r2{hM) 

rc Ir^li 



h ki I2 

. j: i:(2,, + 1),2,. + 1 2 i > } { ;^ } { t i 

rck2 hU 



{pcrb) 



rck2 



f \X,,{cpar)r^,\h,l2), 

Jr Jp 1-1 J 



a ^ b 



rsk2 I3I4 



xXu,{pqrs)T:t{hM) + E E(2^i + + 1)(-1) 



cdk2 



ja+Jp+jb+3q+ki+k2+l3+ll j ^2 ja jc 

jp k h 



h h fci \ / k2 jd jb 
h h k2 



(A.5) 
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